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Abstract 

We analyze the collective gravitational interaction among gas clouds in the inner re¬ 
gions of galactic disks and hnd that it leads to accretion at a rate ~ 
where Mmc is the molecular mass of the disk, Mt is sum of the central plus any axisym- 
metrically distributed mass, and Vt is the mean angular speed of clumps. We discuss 
applications of this result to the mega-maser galaxy NGC 4258, for which we have obser¬ 
vational evidence that the maser spots are concentrated in a thin molecular disk which is 
clumpy, and hnd the accretion rate to be ~ 1.5 x 10“^ Mq yr“^. If the gravitational energy 
release of this inward falling gas were to be radiated away efficiently, then the resulting 
luminosity would greatly exceed the observed central luminosity of NGC 4258, indicating 
that most of the thermal energy of the gas is advected with the how into the blackhole as 
proposed by Lasota et al. (1996). 

The gravitational interactions among molecular clouds lying within the inner kpc of 
our galaxy give an accretion rate of ~ 10“^ Mq yr“^, which is consistent with the value 
obtained by Narayan et al. (1995) by htting the spectrum of Sagitarrius A*. We also 
discuss possible application of this work to quasar evolution. 

Subject heading: accretion disks - galaxies: individual (NGC 4258) 


Alfred P. Sloan Fellow 


1 


1. Introduction 


The gas in the inner parsec of a nnmber of active galaxies appears to be clnmpy and 
concentrated in a thin disk. The best evidence for this conies from the observation of 
discrete spots of water maser emission in active galaxies snch as NGC 4258, 1068, 4945 (cf. 
Watson & Wallin 1994, Greenhill et al. 1995, Miyoshi et al. 1995, Greenhill et al. 1996 
& 1997, Braatz et al. 1996, and references therein). The Inminosity of the central sonrce 
in these systems is believed to be dne to the release of gravitational energy of infalling 
gas. The gas accretion conld be a resnlt of tnrbnlent viscosity in the disk, however its 
effectiveness in a cold thin disk consisting of discrete molecular clouds in very uncertain. 

Here we consider another, more likely, process driving gas accretion - the gravitational 
interaction among clumps, or clouds, excites epicyclic oscillation, causing clouds to collide 
(inelastically) and as a result some fraction of them fall inward. We treat clumps as point 
mass when considering their gravitational interactions, which excite epicyclic oscillations. 
The hnite clump size is taken into consideration when determining the frequency of physical 
collisions among clumps. We assume that there is no dipole or low multipole order coherent 
density structure in the disk which induces efficient inward flow of gas. The effect of spiral 
density wave on accretion in NGG 4258 has been considered by Maoz and McKee in a 
recent paper. The process we consider has some similarity to the work of Goldreich and 
Tremaine (1978), however unlike the case of the planetary ring systems considered by these 
authors the gravitational interactions between clouds are very important and is considered 
in some detail here. 

There is another motivation for considering the gravitational interaction among gas 
clumps. The mechanism for the inward transport of gas within the inner kpc of active 
galaxies (AGNs) to supply fuel for the central activity, is poorly understood. Tidal torques 
from neighboring galaxies, even in close encounters, are believed to be ineffective at this 
distance (but it can be important in transporting gas from several kilo-parsecs to within 
about a kpc). Moreover, recent observations of a sample of spiral galaxies show that the 
central activity is uncorrelated with the strength of the galactic bar (e.g. Ho, Filippenko 
and Sargent, 1997), which suggests that bars are not effective in inducing radial inflow of 
gas either. We estimate the efficiency of gas transport, within the inner kpc of galactic 
disks, due to gravitational interaction between clumps in the disk (Shlosman & Begelman 
1987, also investigated this mechanism; see also Paczynski 1978, and Shlosman, Begelman 
& Frank 1990 and references therein). 

We provide an estimate of gas accretion rate in a clumpy disk in §2 and apply it to 
the water maser galaxies and fueling of AGNs in §3. 

2. Accretion in a clumpy disk 

The accretion rate due to gravitational interaction among clumps is derived in §2.2.1. 
Although pressure forces are ignored in this paper, the method of §2.1 can be generalized 
to include pressure of the inter-clump medium. A physical, but crude, derivation of the 
main result is also presented in §2.2. 

§2.1 Basic equations 
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We consider a thin molecnlar disk which contains clnmps of gas. We assnme that the 
disk does not have some nnstable dipole or other low order long wavelength mode which is 
excited to large amplitnde cansing efficient transport of angnlar momentnm ontward and 
enabling gas to fall inward. However, we assnme that some local instability, gravitational 
or thermal, exists and is responsible for the Inmpiness of the disk. The clnmps move 
on nearly circnlar orbits nnder the infinence of an axisymmetric potential '^(r) and the 
gravitational attraction of other clnmps. The eqnation of motion of the clnmps, neglecting 
pressnre forces, is given by 


(fr 

dt'^ 
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where a{t, r) is the snrface mass density of gas, and r is the center of mass of the clnmp. 
Expanding the angnlar dependence of the snrface density cr{t,r, (f)) in a Fonrier series, and 
snbstitnting the spherical harmonic expansion for |ri — r|“^ in the eqnation above, we 
obtain 
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max{r, ri). It can be shown that the width of Km is r/m, and 


dr I rimKm{r, ri) ~ r. 


( 4 ) 


If clnmps were to move on circnlar orbits with angnlar speed 0(r), with no change to 
their density or shape, then the time dependence of am{t, r) wonld be given by exp \imVt{r)t\ 
We explicitly factor ont this time dependence and nsing eqnation (2) find the torqne, T, 
on a clond of mass to be 


T = 


OO 

27 rmbG m / dri riitrTn(G?'i) exp[tm(/) + . 

__ Jo 
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Note that {4> + Ot) does not change as we follow a blob on a circnlar orbit, however, 
the mntnal gravitational force of clnmps pertnrbs the orbit, and therefore the torqne on 
clnmps finctnates with time cansing the amplitnde of their epicyclic oscillation to nndergo 
a random walk. We estimate the magnitnde of the torqne nsing eqnation (5). 

The torqne on a clond, in a differentially rotating disk, is dominated by nearby clonds. 
To see this, consider two rings of material at radii ri and r 2 , separated by di 2 = |ri — r 2 |. 
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The torque on a cloud in one ring due to clumps in the other ring comes from multipole 
moments am with m < (ri + r2)/2(ii2. The correlation length for am is approximately 
equal to the typical size of a clump, Ij,. Thus the contribution from low order multipoles 
to the torque is small because within the width of the function Km, rijm, we have a 
large number of uncorrelated rings which give positive and negative torques that tend to 
cancel each other outd Furthermore, the correlation time for torque in a shearing disk 
falls off with distance as 0 “^((ib/di 2 ) |(ilnO/(ilnr| (where df, is the average separation 
between clouds and Q{r) is the angular velocity). So clouds lying at radial separation 
much greater than dj, are very inefficient in transferring angular momentum, and most of 
the contribution to T comes from nearby clouds. 

§2.2 Estimate of velocity dispersion and accretion rate 

The number of clouds in a ring of radius r and thickness is ~ (r/b/d^), and therefore 
o'm ~ , where at is the mean surface mass density of clumps and It, is their 

mean size. It thus follows from equation (5) that the torque on a cloud is 


dt) 

where d{r) = abilb/dhY is the mean surface mass density of gas at r, and rumax ~ t/ db- 
The correlation time for the torque tcor ~ ^~^{db/dt^] where dt ~ r{mb/MtYf^ is the 
tidal radius of the cloud, and Mt is the total axisymmetrically distributed mass contained 
inside the radius r. Thus the change to the angular momentum of a clump in time tcor is 
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and the corresponding changes to the clump’s radial excursion amplitude and the radial 
velocity are given by 
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where = Trr^d is the total cloud mass inside of radius r. We expect clouds to undergo 
collision on time scale of tcor when G is not much smaller than dt, and in this case the 
cloud velocity dispersion is equal to dVr given above. The expression for dVr, in the above 
equation, is same as that given by Gammie et al. (1991) for the velocity dispersion of 
molecular clouds in the Galaxy. 


^ We are considering the case where there is no coherent large scale dipole or other low 
order anisotropy in the disk density distribution. 
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Equations (8) & (9) can also be obtained by considering the gravitational interaction 
between clumps and the resulting epicyclic oscillation. The change to the displacement 
amplitude of the epicyclic oscillation of a clump as a result of gravitational interaction with 
another clump, with impact parameter ofd> dt, can be shown to be ~ rubr^/[Mtd?]. Here 
the impact parameter is the distance of closest approach for the unperturbed orbits of the 
two clumps. Gravitational encounters with d < dt, are almost adiabatic and there is little 
change to the epicyclic energy of clumps (two particles moving on circular orbits merely 
interchange their trajectories in such an encounter). Thus the dominant gravitational 
interactions for exciting epicyclic oscillation are those with impact parameter of ~ dt, and 
the change to the epicyclic amplitude in such an encounter is 5r ~ mbr^/{Mtdj) ~ dt, 
which is same as the expression in (8). These encounters occur on a time interval of 
Vt~^{db/dtY ~ icon cind so long as the cloud size {lb) is not much smaller than db the 
collision time is also of the same order i.e. clouds undergo physical collision on the average 
after having undergone one strong gravitational encounter, and so the velocity dispersion 
is ~ Odr, which is same as equation (9). 

We have neglected the drag force of the interclump medium in our calculation which 
is justihed because the turbulent ram force is small compared to the gravitational force so 
long as the density of the inter-clump medium is less than ct. Or in other words the drag 
force is small when the density contrast between clump and the interclump medium is a 
factor of two or greater. 

The cloud collisions are highly dissipative. We assume that the kinetic energy of 
relative motion of colliding clouds is dissipated and their orbit is circularized. Moreover, 
colliding clumps coalesce so long as their size is not so large that rotation prevents the 
merger. Clouds are likely to fragment when their size exceed the maximum length for 
gravitational instability, Imax ~ /Mt, preventing them from continuing to grow in 

size when they undergo collision. In any case there is empirical evidence that gas in our 
galaxy, and in the disk of mega-maser galaxies, is clumpy (see §3). The case considered 
below is when collision between clumps are frequent so that the amplitudes of their epicyclic 
excursion do not grow to be too large. 

§2.2.1 Accretion rate 

The characteristic time for a clump to fall to the center can be estimated using equation 
(8) and is given by 



and the average mass accretion rate is 


( 10 ) 


As pointed out earlier, this result applies so long as the density contrast between the clump 
and the inter-clump medium is about a factor of two or larger, and h is of order db- 
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For Ifj db, clouds undergo several strong gravitational encounters with other clouds 
with impact parameter ~ dt, before undergoing a physical collision, and in this time dura¬ 
tion, t^Qll Vt their epicyclic velocity amplitude, or the velocity dispersion, 

becomes 


5Vr ~ bidt 



(rO) 


db 


(r4)i/2 



and the resulting mass accretion rate is 


M ~ {ftMr) 


db 



( 12 ) 


(13) 


Since the collision time in this case is much greater than 0“^, the effective viscosity is 
suppressed by a factor of {VttcoiiY (^^e Goldreich and Tremaine 1978, for a discussion). 
The physical derivation given above automatically included this factor. Equations (12) & 
(13) are inapplicable in the limit lb 0, when the collision time becomes very long and 
epicyclic amplitudes saturate at a value so that the Toore Q-parameter is about 1. 

For /fe ~ db ~ Imax, the accretion rate in both of the cases considered above reduces 
to 


2 

and the cloud velocity dispersion is given by 

v^rQ . 

\7r MtJ 

Note that the relative velocity of collision between clouds is smaller than the orbital speed 
by a factor of the ratio of the total mass to the molecular mass, and for Mr Mt 
collision speeds are marginally supersonic. The cooling time for gas depends on the density, 
composition and temperature, and it must be short or comparable to the time between 
collisions in order for the system to be in steady state; a specihc case is discussed in §3.2. 

We see from equation (10) that for db ~ Imax the effective kinematic viscosity is 
approximately l^ax^ ~ where Q ~ Imax/Hz is the Toomre Q-parameter, Hz ~ 

v/Vl. This effective viscosity is same as in the ansatz suggested by Lin & Pringle (1987) 
for self-gravitating-disk. We also note that the effective Q of the disk, for the velocity 
dispersion given in equation (15), is of order unity, and thus our solution is internally self 
consistent. 


(14) 

(15) 


3. Some Applications 

In this section we apply the results of §2 to molecular clouds within the inner kpc 
of our galaxy (§3.1), molecular disks in mega-maser galaxies (§3.2), and the evolution 
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of AGNs (§3.3). In all these cases, we estimate an accretion rate and compare it with 
previonsly proposed specihc model of the system. 

§3.1 Molecnlar clonds in the Galaxy 

As a hnal example we consider molecnlar clonds within the inner kilo-parsec of onr 
own galaxy. The amonnt of gas within the inner kiloparsec of onr galactic center is ~ 10® 
Mq. From the molecnlar and HI line emissions we know that the gas distribntion is very 
Inmpy and mnch of it is contained in molecnlar clonds of radins abont 10 pc. The total 
stellar mass within the inner kiloparsec is abont 2x10^® Mq. Thns nsing eqnation (14) 
we estimate the accretion rate to be ~ 10“® Mq yr“^ which is consistent with the valne 
of 10“^a/year obtained by Narayan et al. (1995) from a ht to the observed spectrnm of 
Sagitarrins A*, from radio to hard x-ray wavelengths, in their advection dominated flow 
model. 


§3.2 Mega-maser galaxies 

Water maser emission has been detected from abont a dozen galaxies, of which NGG 
4258 is perhaps the best stndied case (Braatz, Wilson & Henkel, 1996). The maser spots 
in NGG 4258 fall almost along a straight line in the sky from which it is inferred that these 
sonrces he in a thin disk the width of which is observationally estimated to be less than 
abont^ 0.003 pc, and the radii of the inner and the onter edges of the disk are 0.13pc and 
0.24pc respectively. From the observed acceleration of masing spots the central mass is 
inferred to be 3.6x10^ Mq (Mayoshi et al. 1995, Greenhill et al. 1996). The dynamical 
range of observation for NGG 4258 is abont 10® and the maser amplihcation factor is 
estimated to be ~ 10^ (cf. Greenhill et al. 1995, Herrnstein et al. 1997), from which 
we infer that the density contrast between the maser spots and the mean disk is abont a 
factor of three.® The linear drift in the observed velocities of the systemic spots in NGG 
4258, over a period of almost 10 years, is direct evidence that the masing spots are discrete 
clnmps which maintain their identity as they move on circnlar orbits (Haschick et al. 1994, 
Greenhill et al. 1995, Miyoshi et al. 1995). The density of the clonds mnst be at least 
1.5x10^^ cm“® in order that the masing spots are not disrnpted by the tidal held of the 
central mass. This is less than the maximnm density of 10^^ cm“® at which the npper 
energy state of H 2 O is collisionally de-excited and the maser amplihcation is qnenched (cf. 
Reid & Moran 1988, Greenhill et al. 1995). The total molecnlar mass of the disk, nsing 

^ The disk thickness is expected to be ~ 0.002 pc if the molecnlar gas at a temperatnre 
of 10® K is in vertical hydrostatic eqnilibrinm. 

® There is considerable nncertainty in estimating the density contrast from maser emis¬ 
sion becanse the amplihcation factor for the extragalactic masers, if nnsatnrated, is highly 
nncertain (the density contrast wonld be ~ 2 if the gain factor for NGG 4258 were to be 
10®) and becanse the observed maser hnx depends on the path length throngh the disk 
over which the line of sight velocity is nearly constant. A density contrast of 2 to 3 is 
obtained by assnming that the ehective path length does not vary rapidly across a maser 
spot, which is reasonable, bnt for which we do not have independent observational snpport. 
However, we note that if the hnid velocity were to vary on the length scale of a spot size 
then it will necessarily lead to signihcant density inhomogeneity. 
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the above parameters, is thus inferred to be about 2x10^ Mq. The size of the maser spots, 
~ 10“^ pc, is approximately the maximum allowed size for gravitational instability (Imax), 
which is consistent with the suggestion made in the last section that clumps grow in size 
when they undergo collision until fragmentation limits their size to ~ Imax- 

Using equation (15), we hnd that the mean random velocity of clumps in NGC 4258 
molecular disk should be ~ 2 km s“^, which is consistent with the scatter in the ob¬ 
served velocity after the best htting Keplerian velocity has been subtracted from the data 
(Mayoshi et al. 1995). Cloud collisions at this speed are barely supersonic (the gas tem¬ 
perature is expected to be between 500 and 10^ K) and molecules are not dissociated, nor 
do we expect much density enhancement in such a collision. The time between collisions 
can be shown to be about times an orbital period provided that It, ^ Imax- 

Thus we expect a cloud to undergo a collision once every few orbits; these collisions last 
for about Vt~^{lbMt/rMr) ~ a few orbits (here lb is the cloud size). Colliding clouds lose 
energy via rotational and vibration transitions of several molecules of which H 2 O and CO 
are particularly important for the NGC 4258 disk. Using the cooling functions computed 
by Neufeld and Kaufman (1993) we hnd that the H 2 O rotation transitions dominate the 
cooling at temperature of ~ lO^K and H 2 number density of ~ 10^^ cm“^, inspite of the 
fact that these lines are optically thick, and the resulting cooling time for colliding clouds 
is a few tens of years which is smaller than the time between collisions. The average energy 
loss rate resulting from cloud collision in NGC 4258 disk is ~ 10^® erg s“^. 

We estimate the mass accretion rate for NGC 4258 using equation (14) and hnd it to 
be 1.5x10“^ Mq per year. If the entire gravitational potential energy release associated 
with this mass accretion rate is radiated away then the resulting bolometric luminosity 
should be about 2x10^^ erg s“^, which is much greater than the observed x-ray luminosity 
in 2-10 kev band of 2x10^^ erg s“^ and the upper limit to the optical luminosity of 10^^ 
erg s“^. However, the observed luminosity is consistent with the accretion rate estimated 
here if most of the gravitational energy release of the infalling gas is advected into the 
black hole as suggested by Lasota et al. (1996). For a general discussion of the advection 
dominated accretion how model, see Narayan and Yi (1994 & 1995). 

The H 2 O maser emission from NGC 1068 is conhned to a disk between radii 0.65 pc 
and 1.1 pc, and the mass of the central blackhole is estimated to be about 1.5x10^ Mq 
(Greenhill & Gwinn, 1997); the uncertainty in these parameters, however, is larger than 
for the NGC 4258 system. If the density of the molecular gas in the disk is similar to that 
of NGC 4258, then the disk mass should be about 10® Mq and the peculiar velocity of 
clumps, determined from equation (15), comes out to be ~ 10 km s“^ which is close to the 
scatter observed in the velocity of the maser knots. Equation (14) yields an accretion rate 
of 0.05 Mq yr“^, which is consistent with the observed total luminosity of ~ 4 x 10^^ erg 
s“^ (Pier et al. 1994) provided that the gravitational energy release of the gas is efficiently 
radiated away. 

§3.3 Activity in galactic centers 

The nuclear activity in the centers of galaxies (AGNs) is believed to be caused by gas 
accretion onto a massive central black hole. There is some evidence that star formation 
and nuclear activity in the inner regions of galaxies can be triggered or enhanced by 



the presence of a close companion. Close tidal interactions can remove snfficient angnlar 
momentnm and enable gas to fall from the onter parts of the galaxy to within abont 1 kpc 
of the center. However, tidal interactions are ineffective at smaller distances. If the gas 
within this region is cold and clnmpy, as for instance is the case in onr own galaxy and 
some qnasar hosts^, then, as discnssed in §2, gravitational interaction among clnmps canse 
some fraction of the gas to fall inward. Since the accretion rate dne to this process scales as 
QMr{Mr/Mt)‘^ (where O is the angnlar velocity at ~ kpc, and M^. & Mt are the molecnlar 
and the total mass contained within a kpc), the accretion decreases rapidly as the gas is 
depleted provided that the standard tnrbnlent viscosity is nnimportant at ~ kpc scale. 
As the fractional gas mass decreases with time (part of the gas falls into the blackhole 
and some fraction is converted into stars), the accretion rate falls signihcantly below the 
Eddington valne and in the advection dominated regime where most of the gravitational 
energy release is carried with the gas inside the horizon (Narayan and Yi, 1995) cansing 
the Inminosity of AGNs to evolve more rapidly than the cnbic dependence of accretion 
rate on disk mass. Assnming that the star formation efficiency in the central kpc of active 
galaxies is same as the giant molecnlar clonds in onr galaxy, we End the time scale for 
qnasar Inminosity evolntion to be abont 5x10® yrs. 

The rapid evolntion of the AGN Inminosity fnnction with redshift conld perhaps be a 
conseqnence of the accretion dominated by the gravitational interaction between clnmps. 
Since the accretion rate, dne to this process, falls off as inverse sqnare of the blackhole mass 
we expect qnasars with more massive black holes to have a shorter activity life, provided 
that gas is converted into stars at a fixed rate and therefore more of the gas mass is nsed 
np in stars instead of fneling the central activity. 

4. Conclusion 

We have calculated gas accretion, in a cold clnmpy molecnlar disk, dne to gravita¬ 
tional interaction among the clnmps which causes fluctuating torque on clouds. Molecular 
clumps, unlike stellar systems, undergo frequent highly inelastic collisions which damps 
their epicyclic oscillation and causes a net inward flow of gas. We find that the accre¬ 
tion rate due to gravitational interaction among clumps is VLMmc{MmchMtY when 
clumps have size ~ Imax] here Imax is the maximum length for gravitational instability in 
a shearing disk, M^c is the gas mass of the disk and Mt is the sum of the central and disk 
mass. 

We have applied this result to water mega-maser galaxies, and hnd that the accretion 
rate for NGG 4258 should be about 1.5x10“® Mq yr“^ which supports the advection 
dominated accretion flow model for this system suggested by Lasota et al. (1996) i.e. most 
of the gravitational energy release of the accreting gas is advected with the flow instead of 
being radiated away. Maoz and McKee (1997) have arrived at a similar accretion rate in 
their model of NGG 4258 by assuming spiral shock in the circumnuclear disk. 

The total mass in molecular clouds within the central kpc of our galaxy is 10® Mq. 
We hnd that the gravitational interaction among these molecular clouds gives rise to an 

^ Recent GO observation of the cloverleaf quasar at z=2.56 suggests that there is about 
4x10^ Mq of molecular gas within 800 pc of the center (Barvainis 1997). 
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accretion rate of ~ 10“^ Mq yr“^, which is consistent with the valne obtained by Narayan 
et al. (1995) from a ht to the spectrnm of Sagitarrins A* in their advection dominated 
model. 

If the accretion rate in AGNs were to be set by the gravitational viscosity discnssed 
here (at a distance of abont 1 kpc from the center), then their Inminosity wonld nndergo 
vary rapid evolntion, as observed, since the accretion rate is proportional to 

We have assnmed in this paper that the disk is not snbject to large scale dipole or 
other low order harmonic instability that can otherwise efficiently redistribnte angnlar 
momentnm and lead to accretion. Moreover, we have assnmed that clond size does not 
grow monotonically when they nndergo collision; clonds fragment when their size exceeds 
Imax- Thns a combination of merger and fragmentation conld maintain a steady state 
distribntion of molecnlar clnmp size. We have also neglected pressnre forces and other 
hydrodynamical effects of clond interaction. A fnlly self consistent calcnlation, withont 
these assnmptions, is desirable and perhaps will have to be attempted nsing a n um erical 
hydrodynamical simnlation of molecnlar disks. 

I am very gratefnl to John Bahcall, Ramesh Narayan and Charles Gammie for detailed 
comments which led to signihcant improvement in the presentation as well as the scientihc 
content of the paper. I thank Charles Gammie for pointing ont several references inclnding 
the paper of Shlosman and Begelman. 
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